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Abstract 

First-principles density functional calculations are performed to investigate the interplay between 
inplane strains and interface effects in 1x1 PbTiOs/SrTiOa and BaTiOs/SrTiOa superlattices of 
tetragonal symmetry. One particular emphasis of this study is to conduct side-by-side comparisons 
on various ferroelectric properties in short-period superlattices and in constituent bulk materials, 
which turns out to be rather useful in terms of obtaining valuable insight into the different physics 
when ferroelectric bulks form superlattices. The various properties that are studied in this work 
include the equilibrium structure, strain dependence of mixing energy, microscopic ferroelectric 
off-center displacements, macroscopic polarization, piezoelectric coefficients, effective charges, and 
the recently formulated A;_L-dependent polarization dispersion structure. The details of our findings 
are rather lengthy, and are summarized in Sec. IV. 

PACS numbers: 77.84.-s, 77.80.bn, 77.80.-e 
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I. INTRODUCTION 



Ferroelectric (FE) materials have found spread applications in microelectronics such as 
sensors, actuators, transducers, etc. l| In recent years, ferroelectric superlattices have at- 
tracted attention for their promising potential in modifying and tuning the structural and 
polarization properties of FE materials. For instance, when forming superlattices with 
BaTiOa, incipient SrTiOs was found to exhibit strong ferroelectricity. [2] Meanwhile, FE su- 
perlattices grown with desired constituents and/or periodicity provide an important field 
to probe and understand the fundamental physics of ferroelectric materials and related 
properties. |^ Among various FE superlattices, those with ultrashort period are of particular 
interest, since the strong interface effect may lead to some properties in the superlattices that 
are drastically different from those in bulk constituents. In short-period FE superlattices, 
one component significantly influences another, making the material properties interesting 
and less predictable. 

In the study of FE physics, another subfield of importance is to understand the strain 
dependence of FE properties. Inplane strain, caused by either lattice mismatch or external 
stress, alters the interatomic interaction in an anisotropic manner, which often gives rise to 
new physics and/or phenomena. For example, inplane strains have been shown to change the 
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critical temperature of BaTiOa by as large as 500°C. 4| Furthermore, different 
were found to possess very different polarization responses to inplane strain, 
larizations in BaTiOa and PbTiOs were found sensitive to lattice mismatch, 
nevertheless displays a surprisingly weak polarization dependence on the inplane strain. [6 1 
More recently, it was theoretically demonstrated that when FEs are under large strains, the 
Y polarization is to saturate, and this polarization saturation was shown to be a general 
phenomenon applicable for different materials. [7] This finding also leads to a nature expla- 
nation on why polarization in some FEs (not the others) displays a weak strain dependence, 
since the polarization in these FEs is approaching the saturation and thus is less affected 
by the inplane strain. [3] While the strain influences on bulk FEs are amply studied, the 
strain-induced effects in FE superlattices are relatively less understood, however. 

In this paper we intend to address a topic which concerns both of the above two sub- 
fields (namely, FE short-period superlattices as well as strain effects), by studying the 
property changes caused by epitaxial inplane strains in ferroelectric PbTiOs /SrTiOa and 
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BaTiOa/SrTiOa superlattices with short periodic length. The topic is of interest for the 
following reason. In short-period superlattices, the interface plays a far more important role 
than in long-period superlattices and in bulk, and consequently, will significantly alter the 
structural and polarization responses to external inplane strains. The interplay — caused by 
strongly interacting interface and inplane strains — makes the strain effects in short super- 
lattices to differ from their bulk constituents and to be potentially much more complex. The 
above viewpoint emphasizes the differences between FE superlattices and FE bulks. On the 
other hand, FE superlattices must bear some resemblance to the FE bulks, since superlat- 
tices are made of individual bulk constituents. The strain responses of the superlattices thus 
must, to a varied degree, reflect and resemble the properties of bulk constituents. Based on 
these considerations, one key purpose of this work is to conduct a side-by-side comparative 
study of strain-induced effects in short-period FE superlattices and in FE bulks. 

The advantage of a comparative study is rather obvious: by comparing bulks (with no 
interface) and short-period superlattices with strong interface, one is able to obtain a direct 
insight and understanding on the interplay of interface effect and strain effect, and on how 
the existence of one component in superlattices affects the other component under different 
inplane strains. This being said, the comparative study nevertheless is not as straightforward 
as it seems to be, for the following reasons, (i) In certain short-period superlattices, the 
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rotation instability of oxygen octahedral may exist. |8|, 191 On the other hand, bulk PbTiOg 
has only one stable phase of tetragonal symmetry without oxygen rotation as shown by 
the lack of soft modes at the zone boundary [101], and consideration of other phases in 
superlattices makes it difficult to conduct a side-by-side comparison on strain effects between 
superlattices and bulks, (ii) Properties such as ferroelectric off-center displacements, effective 
charges, and polarization structure depend on structural symmetry. By allowing rotation 
instability, most of properties in bulks and in superlattices can not be directly compared, 
and the advantage of comparative study will be largely lost. To enable a comparative study, 
we thus deliberately confine ourself to FE superlattices and bulks of tetragonal symmetry 
without oxygen rotation. For readers who are interested in superlattices with structural 
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phases other than tetragonal symmetry, results can be found in previous reports. [8|, |9|, [111, 



12| Experimentally, superlattices of tetragonal symmetry without oxygen rotation can be 



realized in several possible ways: (1) One can grow short-period superlattice films between 
metallic electrodes possessing no oxygen rotation. The lack of oxygen rotation in electrodes 
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will inhibit the rotation instability in the superlattices. (2) One may use compressive inplane 
strains to suppress the rotation instability. It was shown in SrTiOs that under sufficiently 
large inplane strain, the structure with rotation instability becomes less stable and eventually 
disappears. [3] (3) One may engineer superlattices by choosing atoms of different sizes to 
weaken the rotation instability. For instance, the rotation instability in BaTiOa/SrTiOa is 
considerably weaker than in PbTiOs/SrTiOa. 

In this study, we apply ffist-principles density functional theory (DFT) calculations to 
1x1 PbTiOg/SrTiOg (PT/ST) and BaTiOg/SrTiOg (BT/ST) superlattices, as well as to the 
individual bulk constituents PbTiOg (PT), SrTiOg (ST), and BaTiOg (BT). As shown below, 
various properties are to be investigated, which include microscopic ferroelectric off-center 
displacements, macroscopic polarization, piezoelectric coefficients, effective charges Z*, and 
the recently formulated fc_L-dependent polarization dispersion structure |l3|. A number of 
interesting differences between strain effects in superlattices and those effects in bulks have 
been found, the details of which are summarized in Sec. IV. 

II. THEORETICAL METHOD 

We use first-principles density functional theory within the local density 
approximation (LDA) to determine the structure response to external strains in 
superlattices and in bulks. For the 1x1 superlattices, each unit cell consists of 10 
atoms. The system of tetragonal symmetry has lattice vectors ai, sl2, and sl^, with 
= |a2| = a and |a3| = c. The optimized cell structure and atomic positions are obtained 
by minimizing the total energy. More specifically, for each inplane lattice constant a, the 
out-of-plane c length and atomic positions are optimized. Biaxial in-plane strain is defined 
as rji = rj2 = {a — ao)/ao5 where oq is the equilibrium in-plane lattice constant. In our study, 
we consider compressive inplane strains. 

Calculations are performed using the mixed-basis pseudopotential method. 15|] The norm- 
conserving pseudopotentials are generated according to the Troullier-Martins procedure, [l^ 
Atomic configurations for generating pseudopotentials, pseudo/all-electron matching radii. 



and accuracy checking were given elsewhere ITj. The wave functions of single-particle Kohn 



Sham states in solids are expanded in terms of a basis set which consists of the linear 
combination of numerical atomic orbitals and plane waves. An energy cutoff of 100 Ryd is 
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used throughout all calculations for both bulk materials and superlattices, which we found 
sufficient for convergence. 

In ferroelectric crystals with tetragonal symmetry, the polarization is nonzero along the 
out-of-plane c axis. Polarizations are calculated using the geometric phase of the valence 
manifold of electron wave functions according to the modern theory of polarization [l^ [l^, 
which we have implemented in our mixed-basis computational scheme. 



III. RESULTS 

A. Structure and polarization under zero strain 

To better understand the strain effects in PT/ST and BT/ST superlattices, we first study 
the equilibrium structure of the superlattices under zero external strain. By minimizing the 
total energy of these superlattices, we obtain that the unstrained superlattices have an opti- 
mized inplane lattice constant a=3.87A and tetragonality c/a=2.0 for PT/ST, and a=3.9lA 
and c/a=2.0 for BT/ST. Comparing with our theoretical inplane lattice constants of un- 
strained pure bulk materials — which are a=3. 88A for PT, a=3.95A for BT, and a=3.86A 
for ST, we thus see that, when transitioning from bulk to superlattices, the PT (or BT) 
layers are compressively strained while the ST layers are stretched. 

For unstrained PT/ST and BT/ST superlattices, our calculations using the modern the- 
ory of polarization reveal that both superlattices have zero polarization, showing that the 
properties in the short-period superlattices differ significantly from bulk PT and BT con- 
stituents. The calculated vanishing polarization in the BT/ST superlattice does not contra- 



dict the experimental results |20| where the BT/ST superlattice with one unit cell of BT was 
found to be ferroelectric, since the experimental sample was grown on SrTiOs substrate with 
an inplane lattice constant of a=3.86A. In our case, the BT/ST superlattice is free standing 
without substrate (a=3.9lA). The null polarization in equilibrium PT/ST and BT/ST su- 
perlattices is interesting and meanwhile puzzling, if one recognizes that bulk PT has a very 
large polarization of ~65/iC/cm^ when strained to a=3.87A (the inplane lattice constant of 
the PT/ST superlattice), and bulk BT also has a large polarization of ~35yuC/cm^ when 
strained to a=3.9lA (the inplane lattice constant of the BT/ST superlattice). One may 
wonder why the strained PT or BT layers inside the superlattices do not polarize the ST 
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layers and lead to some polarization. 

To understand why the unstrained superlattices have null polarization, we examine the 
optimized atomic positions at zero strain, which is schematically shown in FigHJ^a). The 
2;- axis is along the superlattice stacking direction, i.e., the direction of the tetragonal c-axis. 
Here we are interested in the relative atomic displacements, rather than the absolute shifts 
of each atom. By placing Pb at the origin, we define a high-symmetry location along the 
z-direction for each atom; more specifically, the high-symmetry locations for Sr and O'^ 
atoms are at z = |, while those of Ti and O2 atoms are respectively located at 2; = | 
and z = ^ along the z-axis. In FigJTJ^b) we show the 2;-axis atomic displacements of the 
LDA-optimized structure with respect to the high-symmetry locations for different atoms. 
We see that the O2 and O2 atoms undergo a notable displacement that is 0.6% of lattice 
length c. Meanwhile, the Ti atom (which shares the same plane as O2) is also displaced off 
the center of high symmetry, but the amount of displacement of Ti differs from that of O2. 
As a result of the different displacements from Ti and O2, our calculations thus reveal that 
there is a local dipole moment for the Ti-02 plane. However, due to the fact that O2 and O2 
atoms move along the opposite directions with equal amount of displacements (so do Ti and 
Ti'), the local dipole moments of the Ti-02 plane and Ti'-Og plane thus cancel, leading to a 
vanishing total polarization. Since both O2 and O2 shift towards Sr, the inversion symmetry 
by the SrO-plane is thus maintained (to the accuracy of numerical certainty). This explains 
why the superlattice has no polarization at zero strain. 

The opposite displacements of O2 and can not be naively attributed to the size dif- 
ference between Sr and Pb, since their numerical atomic sizes do not differ significantly. 
By analyzing the electron density for the optimized structure and for the high-symmetric 
structure, we found that O2 and O2 move towards Sr because of the strong covalent bonding 
between O and Sr. This makes sense since the covalent nature of the Sr-0 bond is stronger 
than that of Pb-0 or Ba-0. It also implies that if one replaces Sr by other A-site atoms 
with less covalent nature, there may be a possibility to produce polarization. Our study 
thus shows that atoms in zero-strain superlattices indeed are displaced off the center, but 
the displacement pattern maintains the plane- inversion symmetry. 
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B. Dependence of mixing energy on inplane strain 

Thermodynamically, when constituents A and B form the A/B superlattice, the mixing 
energy is defined as 

AE{a) = EA/Bia) - Ej,{a) - EM , (1) 

where -Ea/b is the total energy of a 10-atom cell for the 1x1 superlattice, while i?A (or -^^b) 
is the total energy of one unit cell of bulk material A (B). All energies are calculated at the 
same inplane a lattice constant, with the out-of-plane c lattice constant and atomic positions 
optimized for each structure involved in Eq.([T]). Mixing energy tells us the relative stability 
of the superlattice with respect to the individual constituents. A positive mixing energy 
means that the superlattice is thermodynamically less stable and will segregate into pure 
constituents. Of course, this segregation may take a long time due to the possible existence 
of energy barrier. While polarization properties of FE superlattices have been studied, little 
is known thus far about the mixing energy, for instance, (i) What is the value of mixing 
energy for prototypical superlattice such as PT/ST or BT/ST? (ii) How does the change in 
the inplane strain affect the mixing energy? 

Figure [2] shows the mixing energies AE, as well as total energies Ea/b and Ea + Eb, for 
the PT/ST and BT/ST systems at different inplane lattice constants. From Fig|2l we find 
that, (i) For PT/ST, the mixing energy is always positive in the considered strain region, 
revealing that the superlattice is thermodynamically less stable as compared to segregated 
bulk constituents and an extra energy is needed to build the PT/ST superlattice. (ii) 
However, the BT/ST superlattice turns out to have a negative mixing energy at small 
strains, thus being thermodynamically stable energywise. (iii) When the inplane strain is 
small, the mixing energy AE of PT/ST increases in a linear fashion with the deceasing a 
lattice constant. Interestingly, this increase does not last forever; instead AE reaches its 
peak value at a certain inplane lattice and then starts to decline when a is further reduced. 
In fact, we have calculated AE down to a=3.60A (not shown in the figure), which confirms 
the continuous decline of AE. The non-monotonous strain dependence of the mixing energy 
is found true also for BT/ST. (iv) For the strain region considered, the mixing energy ranges 
from 20 to 100 meV per 10-atom cell for PT/ST, while for BT/ST it is within ±20 meV. 
In other words, the mixing energy is large for PT/ST, while being rather close to zero for 
BT/ST. As a result of the small mixing energy, BT and ST are more likely to form FE 
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alloys, which is indeed true in experiments. 
C. Ferroelectric polarization 

Previous studies on the ferroelectric polarization in superlattice largely focused on the 
case that the inplane lattice constant of the superlattice is fixed to be that of substrate 
SrTiOs- Here, our emphasis is slightly different; we examine the polarization in superlattice 
under varied inplane lattice constants, and study how the superlattice responds differently 
(or similarly) to the inplane strain as compared to the bulk constituents. In Fig|3^a) we 
show the total (electronic+ionic) polarizations in PT/ST and BT/ST superlattices, in com- 
parison with the values in bulk PT, BT, and ST. On the one hand, the result in Fig|3]^a) is 
rather trivial — it shows that for a given inplane a lattice constant, polarization in BT/ST 
superlattice is larger than in ST, but smaller than in BT. As a result, polarization in BT/ST 
superlattice never exceeds that in bulk BT of the same a. Similar conclusion is also true 
for PT/ST. On the other hand, a careful examination of the calculation results in FigJHl^a) 
reveals some interesting observations: (1) At ai=3.86A (which is the inplane lattice con- 
stant of typical substrate SrTiOs), BT/ST has a sizable polarization while PT/ST does not, 
although, at this ai lattice constant, bulk PT has a much larger polarization than bulk BT. 
This suggests that, the concept that a stronger FE component (such as PT as compared 
to BT) in superlattice will polarize better the non-FE component (such as ST) does not 
always work. (2) Below and above a2=3.82A, the polarization curve of PT/ST (similarly of 
BT/ST) undergoes an evident change in the slope. More specifically, the polarization rises 
more slowly when a < 02, as compared to the case when a > 02. We recognize that 02 
actually coincides with the critical inplane lattice constant where bulk ST starts to become 
ferroelectric. The critical property of hulk ST (namely, becoming FE at 02) is thus also 
reflected in superlattices. As the ST component turns ferroelectric when a < a2, PT/ST or 
BT/ST becomes less incipient, thus showing a smaller strain-induced polarization enhance- 
ment. (3) When a > 02, polarization in PT/ST is smaller than in BT/ST. However, when 
a < 02, a crossover occurs, and polarization in PT/ST becomes larger than in BT/ST. 

When bulk constituents A and B form an A/B superlattice, one can use the polarizations 
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of bulk FEs and define an average of polarization 

for a given inplane a lattice constant, where Pa and ^a are the polarization and cell volume 
of the A constituent, respectively. This definition is based on the fact that polarization 
itself (namely, dipole moment per unit volume) is not an additive thermodynamic quantity 
and must be weighted by volume. We then compare the polarization Pa/b of the super- 
lattice (calculated using optimized structure and modern theory of polarization) with P by 
examining AP = Pa/b — P- P can be viewed as the anticipated polarization when one 
combines bulk A and B constituents together into a heterostructure, each with the same 
inplane lattice constant a, but without interaction between them. The AP quantity thus 
reflects mainly the interfacial effect on the polarization, caused by various interactions such 
as the polarizing (or depolarizing) field and size effect. 

The P and AP quantities are given in Figl^b) for PT/ST and in Figj3]^c) for BT/ST. 
For PT/ST, we see in FigjSj^b) that (1) when a decreases, P of the superlattice increases 
faster than P; (2) AP is negative in the strain range considered, namely the polarization 
of superlattice does not exceed P; (3) At small compressive strains, Pa/b and P differ 
significantly. But at large compressive strains, they become close. For BT/ST in FigJ^c), 
AP is negative at small strains. However, when a < 3.84A, AP becomes positive, revealing 
that P of short-period superlattice can in fact exceed the average P of bulk constituents. If 
we examine the magnitude of AP (a quantity that indicates the gain of polarization when 
two materials form superlattice), we see that for PT/ST, the gain runs from —40 fiC/cm'^ 
to zero, while for BT/ST, the gain stays within ±10 fiC/cm'^. 

Generally one tends to think that polarization in superlattice is to be enhanced with 
respect to the average of single materials. This need be taken with caution. As shown in 
FigOl^b) for PT/ST, the polarization of the superlattice is considerably smaller than the 
average P of the corresponding single materials at the same inplane a; while for BT/ST, the 
gain of polarization varies from negative to positive as strain increases. As an outcome of the 
competing effect between polarizing field and different covalent strengths of different A-site 
atoms, the gain of total polarization in the superlattice is thus a collective result influenced 
by the properties of single materials, their interaction, and external strain. 

Technologically, one possible advantage of forming FE superlattices is to tune material 
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properties. Here we examine how piezoelectric coefficients may be tuned in PT/ST as 
compared to bulk PT. We are interested in the proper piezoelectric coefficient 631 = — ^^c, 
where x is related to the c-axis polarization by P3 = fiXC- 631 reflects what magnitude of 
polarization enhancement can be achieved by applying the inplane strain rji. In Fig|3]^d) we 
show the X polarizations in PT/ST and in bulk PT. Fitting the x values over the considered 
strain range yields piezoelectric coefficient 631, and the resulting 631 values are given near 
the fitting lines in FigJSt^d). One sees that, the piezoelectric 631 = 19.1 C/m^ coefficient in 
PT/ST is much larger than the value 631 = 10.6 C/m^ in bulk PT. 



D. Microscopic insight: atomic displacements 

To understand microscopically how different layers in superlattices interact and how they 
collectively respond to inplane strains, we report and analyze in this section the atomic 
displacements occurring in different layers of the superlattices. 

Bulk ferroelectric perovskite ABO3 of tetragonal symmetry consists of two layers — the 
AO layer and the BO2 layer — alternating along the c-axis. The total electric polarization of 
the solid could be viewed as the local dipole contributions from the two individual layers, 
as demonstrated by the Wannier functions and local polarizations in Refj! 
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The relative 



displacements of the A and B atoms with respect to the oxygen centers of the same layer 
are thus important quantities that reveal the origin and amplitude of the polarization in 
the material. In 1x1 PT/ST superlattice, there are four non-equivalent atomic layers, 
namely, the Pb-Oi, Ti-02, Sr-O'^^ and Ti'-02 layers (see FigH^). For the convenience of 
discussion, we define the 2;-direction relative displacement of the cation with respect to 
that of oxygen in each layer as ^z{PhOi) = z{Ph) — z{Oi), Az(Ti02) = z{Ti) — ^(02), 
Az{SrO[) = z{Sr) - z{0[), and Az{Ti'0'^) = z{Ti') - z(O^), where z{A) is the ;z-axis 
position of atom A. Az's in the theoretically optimized structures of PT/ST and BT/ST 
are shown in FigJU where the corresponding displacements in pure bulk materials are also 
given for comparison. 

Let us look at PT/ST ffist. It is known in bulk PbTi03 that Pb has a considerable off- 
center displacement. As a result, PbTiOs is a rather strong A-site FE. In comparison, bulk 
SrTi03 has less ferroelectricity from the A-site. This is indeed con&med by our calculation 
results of /\z{AO) for bulk PT and ST (see the dotted lines in FigJl^). However, in PT/ST 
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superlattice, the displacements are remarkably close for the Pb-Oi layer and for the Sr- 

layer (see two solid lines in FigJl^). Also note that ^z{SrO[) in PT/ST superlattice is 
much larger than the counterpart in bulk SrTiOs, for a fixed inplane lattice constant. These 
results are interesting and tell us that, by forming a superlattice, the SrTiOa component 
becomes a much stronger A-site FE, as compared to bulk ST. Regarding Az{PbOi) [or 
similarly Az{SrO[)], we further recognize that this quantity should be identical to zero if 
the 1x1 superlattice has a mirror inversion symmetry by a plane perpendicular to the 
c-axis. On the other hand, once the inversion symmetry is broken by the appearance of 
ferroelectricity, Az{PbOi) becomes nonzero. This is indeed verified by our numerical results 
in FigJD^a), where Az{PbOi) is zero for a > 3.86A and nonzero for a < 3.86A. We thus 
see that Az{PbOi) serves as a microscopic order parameter for ferroelectricity in the 1x1 
superlattice. And this microscopic order parameter can be probed using x-ray diffraction 
since it is atomic displacement rather than electrical polarization. 

In FigJD^b) we examine the Ti relative displacements, Az{Ti02) and Az{Ti'0'2), with 
respect to oxygen in PT/ST. Unlike the A cations where Az{PbOi) and Az{SrO'i) are 
almost identical in different layers, the Ti-0 relative displacements in two Ti02 layers are 
evidently different in FigJDo. At zero strain (a=3.87A), since the O2 atom moves up and the 
O2 atom moves down due to the fact that the Sr-0 bond has a stronger covalent nature than 
the Pb-0 bond as described in a previous section (see FigJT^), Az{Ti02) and Az{Ti'0'^ 
appear to be equal but with opposite sign. With increasing strain, the O2 atom starts to move 
downwards as ferroelectricity is developed, which causes Az{Ti02) to change from negative 
to positive in FigJl^. Interestingly, even for very large inplane strains (e.g., at a=3.75A), the 
difference between Az{Ti02) and Az{Ti'0'2) still exists, showing that the stronger covalent 
nature of Sr-0 bond continues to manifest itself in the microscopic picture. 

From FiglD^a) and (b), one thus sees that, even at large compressive inplane strains, 
the relative atomic displacements in PT/ST are considerably smaller than the counterparts 
in bulk PT, and meanwhile much larger than in bulk ST. This demonstrates the strong 
influence between two constituents when they form superlattice. On the other hand, within 
the PT/ST superlattice, atomic displacements are rather uniform in different layers, except 
for the slight difference in Ti-0 displacements caused by the different covalency in A sites. 

We next examine the situation in BT/ST as shown in FigJU^c) and (d). At small strain 
in FigJD^c), the relative displacements of Ba-Oi and Sr-O'^^ are small and similar. As strain 
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increases, the difference increases. This is in difference from what we have previously seen 
in PT/ST where Az{PbOi) and Az{SrO'i) are close over a wide range of considered strains. 
The difference could be attributed to the large size of Ba atom which, as the inplane a 
constant decreases, will push away more strongly the oxygen atom on the BaO plane, leading 
to a larger difference in Az{BaOi) than in Az{SrO[). Regarding the Ti-0 displacement in 
BT/ST, we see in FigJU^d) that, as compressive strain increases, Az(Ti02) and Az{Ti'0'^ 
become gradually close to each other. At a=3.75A, Az{Ti02) exceeds Az{Ti' 0'2) ■, showing 
a crossover that does not occur in PT/ST. 

E. Effective charges 

In this part of the section we study how the effective charges of atoms are modified when 
forming superlattices. For each inplane lattice constant, we compute the effective Z33 by 
finite difference Z33 = where Ar^ is chosen to be 0.002c. All effective charges are 

given in unit of one electron charge. 

For equilibrium structures of zero strain, the calculated effective charges in superlattices 
and in bulk materials (each at its own equilibrium) are given in Table [B The most notable 
results in Table [T] are: (i) ZI3 of Ti atom in bulk PT is merely 5.76. However, its value 
drastically increases to 7.31 in the PT/ST superlattice. Therefore, Z*s in bulk materials 
can not and should not been used in superlattices. In bulk PT, the Ti atom is strongly 
bonded to only one of the nearby Oi atoms due to the strong tetragonality. In PT/ST, the 
Ti atom is bonded to both Oi and O'^^, leading to a large Z^^. (ii) while Z^^s of the A, Ti, 
or Oi sites change significantly from bulk to superlattice, Z33S of the O2 site are similar in 
bulk and in superlattice. (iii) In PT/ST at equilibrium, Z*s of two non-equivalent Ti atoms, 
namely Zl^{Ti) and Z^^iTi'), are identical, so are ^33(02) and ^33(02). This is caused 
by the planar inversion symmetry in equilibrium structure. Similar conclusion is true for 
BT/ST except for some numerical uncertainty, (iv) In PT/ST, ZI3S of two non-equivalent 
Oi atoms — i.e., ^33(01) and ^33(0^) — are very different. The O'^^ atom on the SrO layer has 
a much larger Z33 than the Oi atom on the PbO layer. 

Under the application of strain, effective charge for each atom is given in FiglSt^a) for 
PT/ST. First, we see that, as the inplane lattice constant decreases by 0.14A, the effective 
charge of Ti atom decreases sharply from 7.31 at a=3.88A to ~4.50 at a=3.74A, demon- 
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strating a wide range of tunability. Similar scale of tunability also occurs to the Oi and O'l 
atoms. In contrast, the Z33 charges of Pb, Sr, O2 and O2 atoms subject to relatively small 
changes (around 0.3). Second, under compressive inplane strains, the Z33 charges of the O2 
and O2 atoms are no longer identical, unlike the zero-strain case. At zero strain, the O2 and 
O2 atoms are symmetric due to the planar inversion symmetry. Under compressive strains, 
the symmetry between O2 and O2 is broken as ferroelectricity develops. Meanwhile, the 
Ti-Oi bond is considerably weakened as a result of the increasing tetragonality, which is re- 
sponsible for the sharp dechne of Z^^(Ti). Third, we recognize that the absolute magnitude 
of the ZI3 charge decreases for most of atoms such as Pb, Ti, Oi and O2, once the impressive 
strain is applied. One exception is Sr. As a decreases, Z^^ of Sr increases instead, probably 
due to the strain strengthened Sr-Og bond. Many of these conclusions are also true for the 
effective charges in BT/ST which are shown in Fig|5](b), except for two evident differences: 
(1) As a decreases, Z33 of Ba increases, unlike Pb; (2) ^13(02) and -^33(02) in BT/ST are 
very similar over the considered strain range. 



F. Polarization structure 



Polarization structure 13|] reveals how the geometrical phase of individual kj 



string contributes to 



_;he electronic polarization Pei, as described by the modern the- 
ory of polarization [isl. Il9 | in the equation Pei = J dkj_(f){k±), where (l){k±) = 
^ Snli /o^" '^^iK'^nfelafjj'I'^nfc)- ^^^^ band structure, the (f){k±) ~ k± polarization structure 
contains various important microscopic insight into the polarization properties. Further- 
more, it was shown that the polarization structure is determined by, and thus can reveal. 



the fundamental interaction among Wannier functions. [13j While the (f){k±) ~ k±_ dispersion 
of bulk ferroelectric has been studied previously, [l3| the polarization structure of FE super- 
lattices remains interesting and unknown. For instance, when bulk BT and ST form BT/ST 
superlattice, the total polarization is known (Figl3K) to decline as compared to bulk BT. 
However, it is not clear at which k±_ points the (pikx) phases suffer more; will the k±_ points 
near the zone center or near the zone boundary suffer most? Also, how is the polarization 
dispersion in superlattice to be affected by the inplane strain? 

The polarization structures of the BT/ST superlattice at two different inplane lattice 
constants, a=3. 86A and a=3.82A, are shown in FigJHl where the counterpart polarization 



13 



structures of bulk BT and ST at the same a length are also made available for comparison. 
Our calculation results in Figj6^ show that, at a=3.86A, the band width of the polarization 
structure in BT/ST is far smaller than that in BT. When transitioning from BT to BT/ST, 
the reduction of the (f){k±) phase occurs mainly near the Xi and X2 points. In other words, 
the (t){k±) phases near the zone boundary are most affected when forming FE superlattices. 

Based on the consideration that (1) bulk BT and bulk ST have very different polarization 
dispersion at a fixed a lattice constant, and (2) the (pikx) phase is inversely proportional to 
the c-lattice length, 1^ one valid approach to compare, at a given k± point, the (j){k±) phases 
in superlattice with those in bulk constituents is to define an average phase as (f){kj_,a) = 
cj^(a)+CB(a) l^J^i^-i-y q)ca(q) + 0B (^±5 «)cb («)] , whcrc Cj(a) is the c-lattice length of bulk i at 
the inplane lattice constant a, and ca/b is the c-lattice length of the superlattice. All 
quantities in the above equation are calculated at the same a lattice constant. is also 
depicted in Figj6l At a=3.86A, we find that in BT/ST can be described rather well 

by 0. However, this is not the case for a=3.82A. In Figj6b, one sees: (i) Though bulk 
ST still has zero polarization with (f){k±)=0 for all k± strings, the 4'{k±) dispersion in the 
BT/ST superlattice is nevertheless notably similar to that of bulk BT. (ii) 4'{k±) in BT/ST 
is considerably larger than the average (f) phase, demonstrating that the strong interaction 
between the BT layer and the ST layer makes the BT/ST superlattice no longer resembling 
the average of two bulk materials. The strong interaction takes place at 02=8. 82A but not 
at ai=3. 86A, since SrTiOs at 02 is near the critical point of becoming ferroelectric, and can 
thus be easily polarized by the electric field arising from the polarization in the BaTiOs 
layer. 



IV. SUMMARY 



Density-functional calculations were performed to study a variety of properties in 1x1 
PbTiOa /SrTiOs and BaTiOs/SrTiOs superlattices of tetragonal symmetry under compres- 
sive inplane strains. An emphasis is placed on the side-by-side comparison of these properties 
in superlattices and in bulks, which is particularly useful in terms of obtaining insight into 
the rather complicated interplay between inplane strains and interface effects. The inves- 
tigated properties include equilibrium structure, strain dependence of mixing energy, ferro- 
electric polarization, microscopic atomic displacements, effective charges, and dispersion of 
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polarization structure. Our main findings are summarized in the following. 

(i) In zero-strain superlattices without oxygen rotation, while atoms are indeed displaced 
off the centers, the displacements nevertheless maintain a plane inversion symmetry. As a 
result, the superlattices show no polarization. The planar inversion symmetry (and thus 
the vanishing polarization) in zero-strain superlattices originates from the strong covalent 
bonding between Sr and O. (ii) The mixing energy is found small and on the order of 20 
meV for the BT/ST superlattice. For PT/ST, this mixing energy is relatively large and 
ranges from 20 to 100 meV in the considered strain region. The small mixing energy in 
BT/ST is consistent with the fact that BT and ST are more likely to form ferroelectric 
alloys, (iii) Under small inplane strains, the mixing energy is revealed to increase linearly 
with the decreasing inplane lattice constant. However, at a certain (large) inplane strain, 
the mixing energy starts to decline with the decreasing inplane a lattice constant (Figj2]). 
As a result, the mixing energy AE exhibits a non-monotonous behavior. 

On ferroelectric polarization under strains, our calculations show: (iv) At the inplane 
lattice constant of SrTiOa substrate (ai=3.86A), BT/ST has a sizeable polarization while 
PT/ST does not, although at this ai lattice constant bulk PT has a much larger polarization 
than bulk BT. This indicates that a stronger FE constituent (such as PT with respect to 
BT) does not always polarize better the non-FE component in short-period superlattices. 
(v) At a2=3.82A, the polarization-vs-a curve undergoes an evident change in slope for both 
PT/ST and BT/ST superlattices, due to the fact that the incipient ST component starts to 
turn ferroelectric. For both superlattices, the polarization rises more slowly when a < 02, 
as compared to the region when a > 02. (vi) The polarization in PT/ST is smaller than in 
BT/ST, when a > 02. But for a < 02, a crossover occurs, and the polarization in PT/ST 
becomes stronger than in BT/ST. (vii) By defining the average polarization P using the 
values of spontaneous polarizations in bulks, we find that AP = Pa/b — P is negative for 
PT/ST in a wide range of considered inplane strains, revealing that the polarization in 1x1 
PT/ST superlattice does not exceed the average polarization P of bulks. On the other 
hand, for BT/ST, AP is found becoming positive when a < 3.84A. AP ranges from -40 to 
fiC/cm'^ for PT/ST, and varies within ±10 fiC/cm'^ for BT/ST. (viii) The piezoelectric 
631 coefficient of PT/ST is calculated to be 19.1 C /iv?, much larger than the value of 10.6 
C/m2 of bulk PT. 

Regarding the atomic displacements in PT/ST superlattice, our study reveals that (ix) 
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in 1x1 superlattice of tetragonal symmetry, Az{PbOi) or Az{SrO[) acts like a microscopic 
order parameter for the appearance of ferroelectricity. Since this order parameter is the 
change in atomic positions, it can thus be probed using x-ray diffraction, (x) The relative 
atomic displacements Az{PbOi) and Az{SrO'^) in PT/ST are found to be very close, over 
a wide range of inplane strains. The large values in our calculated Az{SrO[) indicate that, 
by forming superlattices with PT, the ST component in PT/ST becomes a rather strong 
A-site FE as compared to bulk ST. (xi) The Ti-0 displacements Az{Ti02) and AziTi'O'^) in 
PT/ST differ, however, owing to the fact that O2 and O2 atoms have a tendency to move in 
opposite directions in order to form stronger covalent bonds with the Sr atom. Calculation 
results further show that this tendency continues to manifest itself even at very large inplane 
strains, (xii) For a given inplane lattice constant, the atomic displacements in bulk PT and 
bulk ST are considerably different. However, after PT and ST form a superlattice, the 
displacements in adjacent PT and ST layers are rather uniform, demonstrating the strong 
influence between two constituents, (xiii) In BT/ST superlattice, since Ba and Sr atoms 
have different sizes, Az{BaOi) and Az{SrO'i) deviate significantly from each other at high 
inplane strains. 

On effective charges, the calculation results tell us: (xiv) in PT/ST under zero strain, 
of Ti is 7.31, much larger than Z^^^5.76 in bulk PT. Furthermore, Z^^s of d and 0^ 
in PT/ST are found to be very different, more specifically, Z33(Oi)=-5.36 and -^33(0^)=- 
6.38. (xv) With application of increasing inplane strains, the magnitude of Z33 drastically 
decreases for Ti and Oi atoms, showing a wide range of tunability. Meanwhile, as a decreases, 
1^33 1 of Sr is shown to increase whereas |^33|s of Pb, Ti, Oi and O2 atoms all decrease. 

Finally, the investigation on the polarization dispersion structure demonstrates (xvi) 
when bulks BT and ST form the BT/ST superlattice, it is the (f){k±) phases near the zone 
boundary that are most affected, (xvii) At a=3.86A, the 0(A;j_) phase in BT/ST is revealed 

— * 

to be much smaller than the phase in bulk BT of the same a, and is quantitatively 

close to the averaged phase, (xviii) However, at a=3.82A, the phase in BT/ST 

superlattice is interestingly similar to that of bulk BT, despite the fact that bulk ST still 

shows no polarization at this inplane lattice constant. Furthermore, our calculations show 
that, when bulk ST is near the critical point of becoming ferroelectric, the strong interaction 
between the BaTi03 layer and the SrTi03 layer makes the (l){k±) dispersion in BT/ST no 
longer resembling the average phase. 
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TABLE I: Effective charges Z33 of atoms in PT/ST superlattice (the 2nd column), in BT/ST 
superlattice (the 3rd column), and in bulk PT, BT, and ST (the 4th-6th columns). In 1x1 super- 
lattice, each site has two non-equivalent atoms. Each system is in its own equilibrium under zero 
strain. 



atoms 


PT/ST 


BT/ST 


PT 


BT 


ST 


A site 


3.32 (Pb) 


2.91 (Sr) 


2.75 (Ba) 


2.57 (Sr) 


3.65 


2.79 


2.56 


Ti site 


7.31 (Ti) 


7.31 (Ti') 


7.45 (Ti) 


7.44 (Ti') 


5.76 


7.02 


7.32 


Oi site 


-5.36 (Oi) 


-6.38 (O'l) 


-5.65 (Oi) 


-6.10 (O'l) 


-4.90 


-5.61 


-5.77 


O2 site 


-2.28 (O2) 


-2.28 (O'a) 


-2.22 (O2) 


-2.21 (O'a) 


-2.28 


-2.13 


-2.06 




FIG. 1: (a) Schematic illustration of atomic positions and the direction of atomic displacements 
(by arrows) in the PbTiOs/SrTiOa superlattice at equilibrium. Individual atoms are labeled, for 
the sake of convenience of discussion, (b) Atomic displacements in the LDA-optimized structure 
with respect to the positions of high-symmetry. The displacements are in units of lattice constant 
c. 
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FIG. 2: (a) The total energy E^/b of the superlattice (sohd squares), Ea + Eb of the constituent 
bulks (solid dots), and the mixing energy AE (empty triangles) as a function of the inplane lattice 
constant, for the PT/ST system. E^/^ + Eb are plotted using the left vertical axis, and 

AE is plotted using the right vertical axis. Symbols are the calculation results; lines are guides for 
eyes, (b) The same as (a), but for the BT/ST system. 
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FIG. 3: (a) Total polarizations as a function of the inplane lattice constant, for PT/ST and BT/ST 
superlattices, as well as for bulk PT, BT, and ST. (b) Comparison between total polarization P 
(LDA-calculated) and average polarization P (defined in Eqj2] using the polarizations of bulk 
materials), for the PT/ST system under different inplane a lattice constants. The difference AP = 
P — P is also shown. P and P are plotted using the left vertical axis, and AP is given using 
the right vertical axis, (c) The same as (b), but for the BT/ST system. In (a)-(c), symbols are 
the calculation results, and lines are guides for eyes, (d) The x polarizations as a function of the 
inplane rji strain for superlattice PT/ST, bulk PT, and bulk ST. The numbers given near the 
fitting straight lines are the piezoelectric 631 coefficient. 
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PT/ST BT/ST 




a (A) a (A) 



FIG. 4: Left panel: relative displacements for the PT/ST system; right panel: relative displace- 
ments for the BT/ST system, (a) Relative displacements Az{PbOi) and Az(SrO[) in the PT/ST 
superlattice under different inplane lattice constants (the calculation results are depicted as the 
symbols on the solid lines). The counterpart displacements in bulk PT and in bulk ST under 
different lattice constants are also shown for comparison (by the symbols on the dotted lines), (b) 
Relative displacements Az{Ti02) and Az{Ti' in PT/ST under different inplane lattice con- 
stants (see the symbols on the solid lines). The counterpart displacements in bulk PT and in bulk 
ST are also shown (see the symbols on the dotted lines), (c) Similar as (a), but for the BT/ST 
superlattice. (d) Similar as (b), but for the BT/ST superlattice. 
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FIG. 5: Effective charges Z33 of the non-equivalent atoms in PT/ST superlattice (the left pannel) 
and in BT/ST superlattice (the right panel), as a function of inplane lattice constant. 
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(a) a=3.86A 




X, ^2 \ 2 



FIG. 6: (a) Polarization structures of the BT/ST superlattice, bulk BT, and bulk ST, all at 

the same inplane lattice constant a=3.86A. The average 4> phase (empty circles) is also shown for 
comparison. The kj_ plane of the Brillouin zone is shown as an inset, (b) Similar as (a), but for 
the inplane lattice constant a=3.82A. 



24 



